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Abstract 

We present the formalism of q-stars with local or global U (1) sym- 
metry. The equations we formulate are solved numerically and provide 
the main features of the soliton star. We study its behavior when the 
symmetry is local in contrast to the global case. A general result 
is that the soliton remains stable and does not decay into free par- 
ticles and the electrostatic repulsion preserves it from gravitational 
collapse. We also investigate the case of a q-star with non-minimal 
energy-momentum tensor and find that the soliton is stable even in 
some cases of collapse when the coupling to gravity is absent. 
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1 Introduction 



Boson stars were first discussed in the work of Kaup and Ruffini and Bonaz- 
zola Jj], 0. They investigated the case of a massive, free, complex scalar 
field. Complex scalar fields with a quartic self-interaction have been studied 
in [||, |[]. The case of charged boson stars was thoroughly investigated in 0. 
The general result is that these boson stars are rather small, with radius of 
order of the inverse Compton wavelength. Gravity plays a stabilizing role. 
Another sort of boson stars are the soliton stars. Their feature is that in the 
absence of gravity they remain stable, as simple solitonic solutions to highly 
non- linear Euler-Lagrange equations of scalar fields || [7J. These are very 
large solitons with typical radius ~ Mpjcf)^ and mass ~ M|, z /0q, where O 
is a typical value of the scalar field. These soliton stars resulted from non- 
topological solitons, earlier discussed in M. Another kind of non-topological 
solitons have been studied in the literature called q-balls || . They appear at 

the minimum of the (t//|0| 2 ) ^ quantity in Lagrangians with a global U(l) 
symmetry. In non-topological solitons discussed in || [7j, the fields rotate in 
their internal U(l) symmetry space with an angular velocity u which is very 
small compared to the mass of the free particles. So, the "kinetic" energy 
resulting from the temporal variation of the field can be neglected. This is 
not the case for q-balls or q-stars where this frequency is of the same order as 
the mass of the free particles, because it is related to the parameters of the 
potential (mass, trilinear parameters etc.) The surface of the non-topological 
soliton stars studied in || [| contains an amount of energy which is of the 
same order as the energy stored in the interior. In contrast, the q-star surface 
contains a negligible amount of energy, because its thickness is of order of 
the inverse Compton wavelength. This condition gives a different relation for 
the total soliton mass and its radius, as we shall see. Q-stars with one and 
two bosons in the global case of the symmetry have been investigated in [TD 



Also, q-stars including a boson and one or more fermions have been studied 
in [[11], |T^, |T^, HI], mainly as a possible descriptive model for neutron stars. 

We expect to observe q-balls in supersymmetric extensions of Standard 
Model |1| [16|, [17| , provided that the scalar potential decreases firstly slower 
than 1 0| 2 and then faster and there is a conserved quantum number con- 
cerning scalar supersymmetric partners of fermions. If this supersymmetric 
theory is an effective, low energy one, approximating a more general, prob- 
ably superstring theory in more than four dimensions, we generally expect 
the scalar field to couple to gravity in a non-trivial way, usually with a \<p\ 2 R 



term, [18, 19], where R is the curvature scalar. Also, in a more realistic the- 
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oretical framework, scalar fields may not be governed by a global symmetry 
but by a local one. So, it would be interesting to investigate the possibility 
of having charged q-stars and to compare their properties when varying the 
field strength or the frequency, which will be proved very important for the 
formation of the field configuration. 



2 The general framework 

We consider a static, spherically symmetric metric: 

ds 2 = e u dt 2 - e x dr 2 - r 2 da 2 - r 2 sin 2 ad(3 2 



(1) 



where v = u(r) and the convention g 00 = e u holds. The action in natural 
units for a general scalar field with a local U(l) symmetry coupled to gravity 
is 



d x\[--g 



R 



16nG 



g^(D^)*(D ul p)-U -\f^ 



where 



where e is the charge (or the field strength). Einstein equations are 



(2) 



K - z 5 » R = 87rGT " 
with the energy-momentum tensor: 

T(j, v = (D^)*(D u <p) + (D^)(D^)* - g^[g^(D aV y(D pV )} 

+9,uU + - A g, v F^F aP - g^F m F vfi . 

The Euler-Lagrange equation for the matter field is 

dU 



(l/(V\9\)D,(V\9\rD,) + -^j <p = 0. 



There is a Noether charge: 



f = y/=9Srxt<p*D ll ip-ipD li ip*) = 



-gg^ v i(ip*dvip - cpd u (p* + 2eA u \(p\ ). 



(3) 



(4) 



(5) 



(6) 
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This charge is conserved according to the equation: 



r» = o. (7) 

The total charge is defined: 



Q = j d 6 xf. (8) 

We want to study static solutions, which means that the metric and the 
energy-momentum tensor must be time-independent, though the matter field, 
ip, may depend on time. A plausible ansatz is to choose 

up = ae lwt . (9) 

This is the ansatz used in order to take q-ball type solutions. Specifically, 
this type of non-topological solitons is stable when a is a step-function and 

1/2 

to = uj crit = (U/p 2 )^. We also choose = (A , 0,0,0) in order not 
to have magnetic fields and in accordance with the static character of the 
configuration. As a simple self-consistency check, we can easily find that the 
energy- momentum tensor, and the metrics consequently, is time-independent. 
With these substitutions the Einstein equations are 

A' = (1 - e A )/r + SnGr 2 e x 
(u + eA ) W + U + a' 2 e- x + (l/2)<e^- A ] , (10) 

v > = ( e x - l)/r + 8nGr 2 e x 
(co + eA ) W -U + <r'V A - (l/2)A' 2 e- u - x ] . (11) 

and the Euler-Lagrange equations for the scalar and the gauge field are 
a" + [2/r + (l/2)(i/ - \')]a' + e x (u + eA ) 2 e~"a - e x ^-a = 0, (12) 

A% + [2/r - (l/2)(i/ + X')]A' - 2ea 2 e x {u + eA ). (13) 

We will now give the formulas for the energy and charge of the configu- 
ration. There are two alternative formulas for the charge: 

POO 

Q = 8n drr 2 {uo + eA Q )a 2 e ix ' u)/2 , (14) 
Jo 
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Q = 47rrM' , r -> oo. (15) 

Also, there are two alternative formulas for the total energy (the mass) of the 
soliton 

POO 

E = 4vr / drr 2 [(u + eA ) 2 a 2 e-" + U + e~ x a' 2 + (l/2)A' 2 e' u - x ]. (16) 
Jo 

If we know the total charge of the soliton we can find the total mass from 
the asymptotic relation: 

e x = (l-2GE/r + GQ 2 /47Tr 2 y\ r -> oo. (17) 



We will mainly use eq. [14], [U| 

We have formulated the general framework concerning the case of local 
U(l) symmetry. In order to examine the global case we put e = and A = 
and the results of |jn| can be obtained. It is also very convenient to define 



6 = uj + eA . (18) 
In the case of the global symmetry 9 — > uj. 

3 The q-star solution 

For the sake of convenience we will use a simple potential of the form 

U = a 2 (l-a 2 + ~<j 4 ). (19) 

We will examine separately three different regions: The interior, the surface 
and the exterior of the soliton. At the interior the metric, gauge and matter 
fields vary very smoothly with respect to the radius, because, as we can see 
from Einstein equations, the metric derivatives are proportional to 87rG times 
ip 4 in rough estimate, which is a very small quantity. Let o"o be a value of 
the order of magnitude of the matter field values inside the soliton. We will 
make the rescallings: 

f = m 0j a = a/a , 6 = 6/a , U = U/a 4 . (20) 

We also define 

B = e-", A = e- X . (21) 
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The idea which the q-ball type soliton rely on, in contrast to other sorts 
of non-topological solitons, is the relation: 



The first part of the condition comes from the relation combining uj cr u to the 

potential parameters (i.e. tOcrit — (£VM 2 )mm)- The second is derived from 
solving the equation of motion of the matter field, so as to find its value. 
Remembering that the gravity becomes important when R ~ GM, where M 
the soliton mass, we can easily find that the star radius is R ~ Mpija\. We 
now define the quantity: 



e = V 8nGa ~ m/Mpi. 



[22) 



This means that e is very small. It is easy to find that the radius of the 
soliton is the same order of magnitude as e -1 , or, to be more definite, if f_ 
is the radius of the interior of the soliton, then for the interior holds 



r = r_x = e kx, 
where k is of order unity. We also define 



< x < 1 



~2 2 ( -1 \ 2 

e = e (e k) . 



(23) 



(24) 



So, for the interior, dropping the tildes and the 0(e) terms the Einstein equa- 
tions are 



l-A-x 



dA 

dx 



2 2 
X K 



2e z \ dx 



(25) 



A 



AdB 
B dx 



2 2 
X K 



9 2 a 2 B 



U-—J^-\ AB 

2e 2 \dx 



(26) 



which result from the Gq = 8tvGTq and G\ = %-nGTl respectively, and the 
Euler-Lagrange equations: 



e 2 B 



dU 

d^ 



= o, 



(27) 



dx 2 



+ 



2 1 /l dA 1 dB\ 
x + 2 \A~dx~ + B~dx~) 



dd 2e 2 a 2 9 
dx A 



(28) 
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These equations are in a manageable form because the derivatives of the 
metrics with respect to x are not too small as those with respect to r and 
the soliton radius in units of do -1 is not "huge" any more. It is a matter of 
simple algebra to verify that eq. ^7] has the correct limit when we study the 
global case (i.e.: 6 — > to) and the gravity is absent (so B — > 1). 

The boundary conditions concerning the soliton interior are A(0) = 1, 
A(l) = 1/5(1) and #'(0) = 0. The first condition results from the freedom 
to rescale the metrics, the second from the fact that outside the soliton (we 
will see that the surface region is extremely thin, so one can define as " matter 
region" the soliton interior) the solution to the Einstein equations should be 
Schwarzschild and the third condition reflects the absence of electric field at 
the center of the soliton. 

Eq. |2?| can provide an exact solution for the matter field, i.e.: 



a 

Then: 



1 + OVB. (29) 



U = -(l+6 3 B 3/2 ). (30) 
3 

The surface is very thin in the case of q-stars and this is a crucial difference 
between them and other non-topological-soliton stars. The surface thickness 
is of order a -1 . In this type of non-topological solitons o ~ ou, in contrast to 
other non-topological soliton. The surface energy is of order a 4 R 2 a~ 1 . We 
saw that R ~ Mpi/a 2 . The energy stored in the soliton interior is Ei nt ~ 
R 3 a A . So the relation: E sur ~ eE int hold. These results are numerically 
verified. So the surface energy is negligible, a result that not hold for the 
other kind of soliton stars. Within the surface region, the metric and the 
gauge field varies very slowly (their derivatives are of 0(e -1 )) and can be 
regarded as constant. The matter field varies rapidly (this is the meaning of 
a surface). Now we will find the "eigenvalue" equation for the gauge field, 
or the frequency in the global case. In eq. [12] 2/r + (1/2) (i/' — A') are 0(e) 



terms and can be neglected. Eq. [12] can be straightforward integrated giving 
as result: 

e- A (J) 2 + ^Ve-- f / = 0. (31) 

In the interior the above equation, neglecting the 0(e) terms, will contain 
the second and third term. In order to confirm the continuity of the above 
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expression (and should be continuous as a first integral) we should set 

(e 2 a 2 e~ u - U) r = 0. (32) 

Switching off both gravity and the gauge field (i.e. making the symmetry 
global) the above expression arises from the Euler-Lagrange equation when 
the spatial derivative terms are neglected (i.e. the soliton is very large) and 
uj = uj min . Now 9 took the place of u. Eq. 152 gives the eigenvalue condition 



for the frequency. Substituting equations |29], |30] into eq. we find that for 
the surface 

9 2 sur B sur = 1/4. (33) 

We can easily understand that different values of the potential parameters 
gives a different relation between 9 and the metrics. When both gravity 
and electric charge are absent the above relation takes the form uj = 1/2, 

which is really the value (U/\(p\ 2 ) min for the potential of eq. 0. Eq. [33| 
gives a measure of the gravity strength when we know the frequency. When 
A(sur) < 1, i.e. when we switch on gravity, the frequency of the large soliton 
(thin-wall approximation) is less than that corresponding to the absence of 
gravity. This fact summarizes the main effect of the gravity upon the soliton 
properties. The "potential" that enters the ui min relation now is not U but 
A(sur)U and is smaller than in the case that the gravity is absent. 

In the exterior region the matter field a is zero. The metric fields should 
obey the asymptotic relation 

A(r) = l/B(r), r — > oo 

One can prove that the above condition holds for the entire exterior region, 
not only at infinity, using equations 

AdB 2 2 1 fd9\ 2 A n 

A-l-x — = -x 2 k 2 — — AB, x > 1 

B dx 2e 2 \dxJ 

which concern the metrics outside the soliton, and the above boundary condi- 
tion. This result is in agreement with the Schwarzschild type of the exterior 
solution. The boundary condition for the gauge field 9 results from the 
matching of the exterior end the interior solutions. The other condition for 
the gauge field is 9 — > uj when r — > oo which means that for large distances 
the "true" gauge field, (9 — u)/e, tends to zero. Our results also hold in the 



global case of the U(l) symmetry studied in (lC . 
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Figure 1: Soliton profiles for different values of the coupling constant, with 
A(sur) = 0.36 or, equivalently, 9 sur = 0.3 for five different values of the field 
strength e. We examine only the interior region because the surface is very 
thin. 

4 Numerical results 

The computation of energy (total mass) and particle number of the field 
configuration are based on eq. [13], [H]. k gives a measure of the soliton radius. 
The main features of the soliton are plotted as functions of the frequency, 
which is the most crucial parameter. Alternatively, we could use a parameter 
of the surface gravity (A or 1/B) in order to test the soliton behavior. But 



9 sur , or u, are not independent parameters, as we can see from equation |33. 
So, the value of the frequency gives a measure of the surface gravity strength 
as well. 

The total particle number can be identified with the total charge if we 
take as unity the charge of each particle. Figure [I] shows the behavior of 
the soliton profile when decreasing slowly the field strength e. We sketched 
only the interior region because the surface is very thin and the field varies 
rapidly. When the coupling constant increases, the soliton gets larger due 
to the electrostatic repulsion and, most interesting, the value at the origin is 
not the maximum one, due to the same effect. 

Fig.|2| shows the behavior of the soliton charge as a function of the fre- 
quency. These results are intuitively expected. When 6 sur or to decreases, 
the particle number increases because the gravity gets more important (eq. 
PB] ). The presence of the electric charge e makes the decrease of the particle 
number slower, due to the electrostatic repulsion between the different parts 
of the soliton. 
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0.1 0.2 0.3 0.4 0.5 

Frequency 



Figure 2: The total charge (particle number) of the soliton as a function of 
9 sur , which equals to uj in the global case, for four different values of the field 
strength. 




Figure 3: The mass of the field configuration as a function of 6 sur for four 
different values of the field strength. 
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0.5 



0.1 0.2 0.3 0.4 0.5 

Frequency 

Figure 4: The soliton radius as a function of 9 sur for four different values of 
the field strength. 

The results of figure ^ are more interesting. In the case of global sym- 
metry, the increase of frequency decreases the total soliton energy due to 
the decrease of the particle number. When the field strength e is large, the 
decrease of the particle number is not enough to decrease the energy and 
the electrostatic energy should be taken into account. This energy becomes 
more important when gravity weakens, i.e. when frequency approaches the 
value that has when gravity is absent. For the potential of eq. ^ this value 

^ (U/lrftl = 0-5. 

The soliton radius, figure f|, has a very interesting behavior. In general 

the radius gets larger with the coupling constant due to the electrostatic 
repulsion between the different parts of the soliton. The radius also increases 
when the frequency decreases, because the particle number gets larger. But in 
the global case and after a certain point the soliton radius decreases because 
the gravity becomes extremely strong. This is the boundary between stars 
and black holes, because when u — > 0, then A(sur) —> (eq. |33] and the 
relation A(sur) = 1/B(sur) between the metrics) and the metric B at the 
center of the soliton approaches infinity. 

We will now study the soliton behavior when varying the coupling con- 
stant e. A general feature is the decrease of the soliton mass, radius and 
particle number when the coupling constant gets larger. We also observe the 
increase in the ratio soliton energy per particle number (which expresses the 
ratio soliton energy per free particles energy) which leads to an "electrostatic" 
instability of the soliton for very large values of the field strength e. The field 



11 




Figure 5: The value of the scalar field at center of the soliton as a function 

of @sur m 




Figure 6: The soliton charge as a function of the coupling constant for four 
different frequencies. 
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gure 7: The total soliton mass as a function of the coupling constant 
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configurations described in the figures satisfy E so i/Ef ree < 1, i.e. they are 
stable configurations. In order to simplify calculations we chose the mass of 
the free particles to be 1 (eq. [19]), so the energy of the free particles with 
the same charge as the soliton equals the particle number. When the above 
ratio exceeds unity, the soliton decays into free particles. On the other hand 
the presence of non-zero electric charge prevents soliton from gravitational 
collapse. In the global case and when uj is very small, the metric B at the 
center of the soliton tends to infinity, revealing a singularity at the center of 
the field configuration. Electrostatic repulsion removes this singularity and 
makes the soliton stable, even for lower frequencies. 

The plot of charge as a function of the field strength, fig. || has a well- 
understood behavior. The increase in the electrostatic repulsion decreases 
the number of particles within the soliton. The plot of energy, fig. |7|, is 
more interesting. When the frequency is small and, consequently the gravity 
is more important (eq. |33[) the increase of the coupling constant leads to a 
decrease of the energy due to the decrease in the total charge. When gravity 
is weaker (for instance 9 = 0.35 or, equivalently A(sur) = 0.49), the increase 
of the coupling constant increases the electrostatic energy. The behavior of 
the soliton radius is interesting as well. Large coupling constant leads to 
bigger solitons when the strengthening of gravity has the opposite effect. 

5 Q-stars with non- minimal energy- 
momentum tensor 

We will now examine the case where the action has a new term concerning 
the coupling of the scalar field to gravity. Such terms can arise in effective 
Lagrangians when considering theories in more than four dimensions. We 
will use a simple, toy model for the action, namely: 



We have regarded only the case of global U(l) symmetry for simplicity. 




(34) 



The energy-momentum tensor is 



T^ u = dp<p*d v <p + d v ip*df,<p - g^(g a Pd a (p*df3(p) + g^U + 
2^>(i^„ - X - 9ilv R) - 2^g^d a dM 2 + ^W<p\ 2 



(35) 
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Using the ansatz of eq. | referring to the static configurations we take 
the Euler-Lagrange equation for the matter field: 



2 1 

- + ^'-A') 
r 2 



z x ~ v u 2 a + e x £Ra - e x ^—a = 0, 



dU 



and the Einstein equations: 



16f7rG) 



A' 



ry ty^ ry 2 



8nG[e x - u LU 2 a 2 + U 



2_2 -1/ 



(36) 



(37) 



v' e x 1 

' ry j>i 



87rG[e A - v a; 2 t7 2 - 



U + a' 2 -2iaa\v' + 4/r)], 



(38) 



[l + 16^7rGa 2 (l + 60] J R = 
8ttG[(2 + 120(e~V 2 - e-"u 2 a 2 ) + 4Z7(1 + 30], 



(39) 



where the last equation has been obtained by taking the trace of the Einstein 
equations and will be used for the calculation of the curvature scalar. 
The total mass of the configuration is 



E = 4vr 



drr z 



, n l + 16£vrG<r 2L V W 

U(l + 40 + e~V 2 (l + 40 - 2e- x £v'a'<j - 4fVi2]. 



(40) 



The total particle number is given by eq. |5j, or, equivalently, [14]. The poten- 
tial is the same as in eq. [L9|. We also make the rescallings and the definitions 



01 and 21. We also define 



R = R/al 



(41) 



We discriminate between three regions: The interior region, the very thin 
surface region and the exterior, where the metric is Schwarzschild. Follow- 
ing the steps of the previous case, we write the equations for the interior, 
dropping the tildes and the 0(e) quantities: 



dA 



I- A- x— = x 2 k 2 \9 2 a 2 B(l - 40 + U(l + 40] , 

dx 



(42) 
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Figure 9: The central field value as a function of the coupling £ of the scalar 
field to gravity for three different frequencies. 



A-l-^ = a*«>[<Po>B-U], (43) 

resulting from the G® and G\ Einstein tensors respectively. The rescalled 
curvature scalar, R is of 0(e) order. So the total mass of the configuration is 

poo 

E = 4tt / drr 2 {e 2 a 2 B{l-AC) + U{l + Ai)). (44) 
Jo 



We also find that eq. [29| for the field value in the soliton interior remains 
valid. 



We solve numerically eq. H2I, E31 and we find the mass and particle number 



of the soliton using eq. H and [14| respectively. Some general features of the 
solutions for the £ ^ cases are: a) For small values of the parameter that 
describes the strength of the coupling of the scalar field to gravity the soliton 
becomes more "stable" against gravitational collapse. We can find solutions 
for some values of the frequency that have B(0) — > oo when £ is zero. These 
solutions describe a real singularity at the soliton origin, which is removed 
when we "switch on" the coupling to gravity. But this coupling, however 
strong, can not prevent soliton from gravitational collapse when A(sur) — ► 0. 
b) Larger values of the coupling parameter correspond to solutions with lower 
energy and charge but larger radius, which means that the coupling of the 
matter field to gravity has an effect similar to the electrostatic repulsion. 
This effect has to do with the repulsive character of the £,<p 2 R term. 



16 



-2 2 4 6 8 10 

Coupling to gravity factor 

Figure 10: The radius of the soliton as a function of the coupling £ for three 
different frequencies. 
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Figure 12: The soliton particle number as a function of £. 
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